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ABSTRACT 


Terminal composition control of a binary distillation column 
operated under computer control has been studied for disturbances in 
feed flow rate. Conventional two point control, whereby the overhead 
composition is controlled by the reflux flow rate and bottoms composi- 
tion by means of reboiler heat duty, was demonstrated to be unsatis- 
factory. The unsatisfactory behaviour of this control system is a 
direct result of the interaction between the two control loops due to 
the dynamic coupling inherent in a distillation process. Two alternate 
control systems, namely a noninteracting control system and a ratio 
control system, whereby the overhead composition is controlled by mani- 
pulating the reflux to overhead vapour rate ratio were proposed in order 
to eliminate the deleterious effects of the dynamic coupling. 

The control systems were evaluated using a 9 inch diameter, 

8 tray pilot scale distillation column interfaced with an IBM 1800 
digital computer. The results show that a very significant improvement 
in the control of both product compositions is achieved by using tne 
noninteracting system or the ratio control system as compared to the 
behaviour that results using conventional two point feedback control. 
However, the control performance obtained using the noninteracting was 
only marginally better than that obtained using the ratio control scheme. 

The noninteracting control system was constructed from a 
simplified transfer function representation of the distillation column 
dynamics, determined from a series of on-line pulse tests. It was 
found that the dynamic behaviour of the distillation column, about an 
operating point, could be adequately represented by a set of first 


order plus dead time transfer functions. 
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1.__ INTRODUCTION 


The distillation process is extensively used in the chemical 
processing industry to separate a feed stream into product streams of 
a specified composition. In order that a distillation column, subject 
to disturbances in a changeable environment, yield the specified 
product compositions, some form of control system is necessary. 

The control of a distillation column to produce a single 
stream at a specified composition has received wide spread application. 
On the other hand, the control of more than one product composition 
presents some special difficulties due to the dynamic nature of the 
process. Consider a two point control system whereby the overhead 
and bottom compositions are maintained at their desired values by 
manipulation of the reflux and steam flows respectively. It is readily 
apparent that manipulation of the reflux flow will affect not only the 
top product composition but also the bottom composition as well. Like- 
wise, the reboiler heat duty will also affect both of the end composi- 
tions. Indeed, the magnitude of the interaction between the two control 
loops, that is the effect of the reflux flow on the bottom composition 
and steam flow on the top composition are of the same magnitude as the 
direct effects of the manipulated inputs on the controlled variables. 
Thus an attempt to control the two product compositions in this manner 
will not meet with success, as control of the top product composition 
will cause disturbances to enter the bottom composition control loop 
and vice versa. 

The synthesis of a control system embodies the appropriate 


selection of the controlled and manipulated variables based on the 
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objectives of the control system as well as the allowable degrees of 
freedom for the process. An analysis of the degrees of freedom for 
a distillation process is presented in Section 1.1-D. If the objec- 
tive of the control system is production of one or more products of 
a specified quality then the main controlled variables would clearly 
be the specified product compositions. However, in many instances 

it is not possible to obtain a direct measure of the product quality. 
As a result, temperatures on intermediate trays, or the column pres- 
sure, or the reflux/feed ratio are often selected as the controlled 
variables. Reflux rate, reboiler heat duty, and product flow rates 


are the most common manipulated variables. 


1.1 LITERATURE REVIEW: 

The control of distillation columns has received much atten- 
tion in the recent chemical engineering literature. It is not possible 
to discuss all the literature pertaining to distillation column control. 
Rather, only those studies relevant to dual quality control, i.e. the 
control of both end compositions will be outlined in detail. For 
descriptions of control systems other than dual quality control, the 
reader is referred to several excellent reviews in the literature (1,2, 


364). 


1.1-A Dual Quality Control 

In a recent series of papers Rijnsdorp (5,6,7) discusses the 
interaction in two variable control systems for distillation processes. 
It is concluded that the interaction inherent in the process will be 
deleterious to controllability in many instances. An alternate control 


scheme, using ratio control between the reflux and top vapour flows is 
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suggested for reducing the effects of the interactions. The top pro- 
duct quality controller is used to adjust the ratio setting. Owing to 
the ratio control top product composition will be insensitive to 
variations in vapour flow, i.e. reboiler heating. This control system 
is discussed in more detail in Chapter 4. 

In a recent theoretical study of distillation column dynamics, 
Davison (74) illustrates that severe interaction may exist in a binary 
distillation control system between temperature, or composition and 
pressure control loops and between the top and bottom temperature or 
composition control loops. 

Distillation column control loops are inherently sluggish due 
to the large time constants within the process. An approach to 
improving the performance of such systems has been suggested by 
Rosenbrock (8) and extended by Gordon-Clark (9), Davison (10) and 
Gould (11). In the formulation presented by these authors, it is 
assumed that the distillation process can be adequately modelled by 


the state equation 

x = AX+BU 
The matrix A may be written 

hes ALL Oly 


where V is the matrix of eigenvectors of A and V' is the matrix of 


eigenvectors of Al. If the control vector u is chosen such that 


u = Ky' x and B = V, then the controlled system becomes 
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Thus the eigenvalues of the system may be made to take any 
desired value by specification of the gain matrix K. If K is a 
diagonal matrix, it can be shown that the shifting of any one eigen- 
value does not influence the other eigenvalues of the system. Thus 
the result will be a form of noninteracting control. 

However, several restrictions will be encountered when an 
attempt is made to apply the above theory. The control matrix B will 
not usually be identical to V. Gould (11) has shown that this is not 
a serious restriction and may be overcome. A more serious restriction 
will be encountered in a system with inaccessible state variables, as 
would be the case in most complex industrial processes. In addition 
in a regulatory control system, it would be necessary to add integral 


control, the integral contributions to control would then interact. 


1.1-B Multivariable Process Control 

The inherent interaction present in a distillation column 
control system would suggest the application of multivariable control 
theory to symthesize a two point control system for the process. There 
are basically two different types of design methods available for 
multivariable control systems. The first method (noninteracting 
control) consists of placing a compensator or feedforward controller in 
the feedback path to compensate for the effects of interaction between 
control loops. The second technique is based on optimal control theory 
and consists of optimizing some objective function given a plant model. 
The effects of interaction are not eliminated in this case rather they 
are taken into account or used to improve system performance by means 


of the optimization process. 
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Synthesis of a Noninteracting Control System 

Dynamic inter-coupling between process variables in a multi- 
variable system prevents the direct application of single variable 
control theory. However, if the interacting effects could be can- 
celled by some form of compensator then single variable techniques 
could be employed to construct a suitable control system. 


Consider the system shown in Figure 1.1. 


FIGURE 1.1 


The closed loop transfer function for this system is given by 


eg ere) 


[real ee By definition the system is noninter- 

acting if the matrix G is diagonal. If G is diagonal, one input will 
effect only one output; hence noninteraction is achieved. It can be 
Shown that if the open loop transfer matrix P C is diagonal then the 
closed loop transfer matrix will also be diagonal. Thus it appears 
that the noninteracting controller can be synthesized by open loop 
methods. The method for open loop synthesis of a noninteracting control 
system as given by several workers (12,13,14) will not be outlined. 


Consider that the elements of the open loop transfer matrix G are given 
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The design procedure involves selecting the elements Cry so that the 
off diagonal elements of G are zero. 
le3 N 
ns Pry Syy = 0 [a7 es 


If the form of the diagonal elements Gry are specified from frequency 
response considerations or other methods, the elements of the com- 
pensator can be computed. The compensator must be chosen so as to 
yield a stable system. Morgan (13) and Chen (14) have outlined a set 
of constraints that must be placed on the compensators to insure 
Stability. In addition to stability considerations the form of the 
Specified open loop transfer function G. must be such that the 
elements of the compensator C are physically realizable. This is a 
trial and error task, indeed the desired form for G. may not be 
possible due to the physical realizability constraint. 

Zalkind (15) gives a method for obtaining noninteracting 
control that is somewhat easier to realize than the above method. 
Tnis method consists of placing two compensators in series as illus- 
trated in Figure 1.2; the compensator D is introduced to cancel tne 
interaction effects and the controller C can then be constructed by 


Single variable methods to control the individual loops. 
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FIGURE 1.2 


In Figure 1.2 C is a diagonal matrix of controller transfer functions. 


The compensator D can be determined from 


DP = diagonal matrix 


and the diagonal elements dij = 1] 


Once again the compensator D may be difficult if not impossible to 
realize depending on the form of P. 

A considerable amount of work has been done on the synthesis 
of multivariable control systems to yield a specified closed loop 
response. The various methods of synthesis when the closed loop 
response is specified have been discussed by several workers (16-28). 
Basically the method consists of specifying a suitable form of the 
closed loop response Y. However, it should be noted that the form of 


Y must be such that the compensator C is physically realizable. 
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solving for C gives 


CRS ae. 


This method has been outlined by Kavanagh (19,20). 

Horowitz (16,22) also presents a synthesis procedure for 
multivariable systems. The objectives of this procedure are rejec- 
tion of disturbances and compensation for changing plant parameters. 

Cther investigators (24,25,26) have discussed the synthesis 
of noninteracting controllers and stability of a multivariable system 
by the use of linear transformations. 

Foster (29) and Greenfield (30) describe methods which may 
be used to decouple a multivariable system by means of a canonical 
transformation. Feedforward and feedback controllers are then added 
to each independent subsystem. 

Kinnen (32) describes a method of multivariable control 
system synthesis based on a root locus analysis. This method is 
particularly applicable to systems where design specifications such 
as rise time, overshoot, etc. are important rather tnan noninter- 
action. Design parameters can be chosen to restrict the interaction. 

As discussed previously, the application of noninteracting 
control theory (11-30) may give rise to an overly complex system. On 
account of this and the fact that it may not always be possible to 
achieve the desired performance, Rosenbrock (33) suggested an alter- 
nate synthesis technique. A structural synthesis technique is out- 
lined which will handle the nonminimum phase effects caused by pro- 


cess dead times and the existence of right-half-plane zeros. 
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Another disadvantage of the previous conventional analysis 
for obtaining noninteracting control is the loss of structural informa- 
tion about the process in the complex form of the controller. Thus a 
design method that would retain the structural characteristics of the 
process may be desirable in some cases. Techniques developed by 
fathias (35), Greenfield (36,37) and Morgan (28) are designed to retain 
Structural information about the process in the controller. 

A multivariable controller synthesis using time domain pro- 
cedures will often be more desirable than a frequency domain approaci. 
This is especially true for a system with a large number of states as 
the frequency domain approach would yield an overly complex controller. 
Morgan (38) and Falb (39) have introduced a technique for attaining 
noninteraction by means of state variable feedback. The method is 
similar in principle to the frequency domain techniques already des- 


cribed. For a process which can be described by 


x(t) = A x(t) +B u(t) 


and if control action of the form 


a saab Cote Salts 
is applied, then the expression for the controlled plant is given by 
Pigs Ube Syne me ental Ns 


Thus for complete noninteraction among the states and the reference 
inputs the matrices [A + B H] and B k must be diagonal. 


Little work has been done on sampled data versions of 
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noninteracting controllers. Fodor (40) has summarized the synthesis 
and limitations of a sampled compensator for noninteracting control. 
The major conclusion of the work is that a system can only be made 
noninteracting at the sampling intervals. 

The synthesis of a decoupling compensator for a distilla- 
tion column is discussed by Luyben (41). The series type compensa- 
tor recommended by Zalkind (15) is suggested for "two-point" control. 
In another article on the multivariable control of phase separation 
processes Kleinpeter (42) illustrates the application of the struc- 
tural formulation of Greenfield and Ward to a vapor liquid separation 
cell. The results serve to illustrate the advantages of a structural 


formulation over the terminal approach for noninteracting control. 


Synthesis of an Interacting Control System 

Since the interaction is an intrinsic property of a multi- 
variable system, some workers have suggested that improved perform- 
ance would result if the interaction is allowed to remain unimpeded 
or used rather than eliminated. 

The majority of methods that appear in the literature 
(43-52) for interacting system synthesis are based on optimal control 
theory. Given a dynamic model of the plant the usual approach is to 
optimize some objective function via dynamic programming or calculus 
of variations. However, it has only been in a more recent study by 
Newell (76) that an optimal control technique has been applied to a 
regulatory control system with unmeasured load disturbances. 


Mesarovic (53,54) has shown that an interacting system can 
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give improved performance over an uncoupled system. From knowledge 
of the disturbance inputs an interacting domain is constructed in 
which the performance of an interacting system is superior to that of 
a decoupled control system. 

Unlike their noninteracting counterparts, optimal control 
systems are more easily formulated in discrete time. Optimal control 
of discretized systems is the subject of much attention in the litera- 
ture. The reader is referred to books by Tou (55) and Lapidus (56). 

In a paper of specific interest Heiderman (52) has illus- 
trated the design of an optimal discrete controller for a system with 
load changes. The technique is applied to a distillation column. 
However, the design method is restricted in that the disturbances 
must be measured. 

A paper by MacFarlane (77), published after completion of 
this study, presents a survey of multivariable feedback control theory. 
This paper reviews the available design techniques and compares the 


advantages and disadvantages of the design methods. 


1.1-€ Process Control Via Gas Chromatograph 

It has long been recognized that a gas chromatograph (GC) is 
a reliable and accurate tool for chemical analysis in the laboratory 
and the same should apply for its use in process control systems. 
Although there is widespread use of the instrument for on stream 
analysis, there is little evidence of process control loops being 


closed with a GC. There are several reasons for this, the discontinuous 
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Signal and measurement lag may seriously affect the controllability of 
the system, however, probably the most serious obstacle is the com- 
putation of the stream compositions from the voltage output of the gas 
chromatograph. 

Villaboo (57,58) outlines the advantages and application of 
gas chromatographs in industry; while Vankampen (21) discusses the 
difficulties of process control via a gas chromatograph and outlines 
the dynamic problems connected with the introduction of a gas chroma- 
tograph analyzer into a process control loop. The effect of sampling 
time on system stability is examined by means of a polar diagram. 
Satisfactory control of a pentane-propylene splitter was obtained with 
a ratio of sampling time to process time constant at .18. Froode 
(59,60) outlines the use of a digital computer - gas chromatograph 
system for process control. The design of the computer system is dis- 
cussed and examples of industrial applications are given. 

The use of a small digital computer to monitor gas chroma- 
tographs surmounts the previously mentioned difficulties of computing 
the compositions and incorporating the results in a control system. 
Dedicated digital computer systems have become popular for monitoring 
laboratory gas chromatographs. Several authors (61,62,63) have des- 
cribed the use of such systems. However, little has been reported on 
an integrated gas chromatograph process control system using a digital 
computer. 

A computer monitor system for gas chromatographs must pro- 
vide for the following functions. It must be possible to initiate a 
gas chromatograph analysis by means of an external or internal inter- 


rupt within the computer. In the case of a process control] system the 
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computer must have the capability of initiating the GC run and control- 
ling the sample injection on a repetitive basis. The system must be 
capable of reading the voltage output from the GC bridge and converting 
the signal to a digital value at rates up to 5 times a second. After 
the data acquisition or portion thereof is complete the monitor system 
must be capable of analyzing the data. Lichtenstein (61) and Briggs 
(62) have outlined the necessary features of the data analysis program. 
The analysis of the data consists of peak recognition, integration of 
the area under the peak and calculation of base line corrected areas. 
The method most commonly used for peak recognition is examina- 
tion of the first and second derivatives. The derivatives may be com- 
puted based on the least squares estimates as outlined by Savitzky (64). 
The start of a peak is recognized when the first derivative becomes 
greater than some specified deadband. The termination of a peak is 
reached after the sign of the first derivative changes from negative 
to positive. Examination of the second derivative is the usual means 
of determining if the base line is reached at the end of a peak. Base 
line corrected areas are computed by saving the coordinates of the 
start and end of the peak and then fitting a least squares line through 
sections of the base line immediately before and after the peak. The 
superfluous areas are then subtracted and the compositions computed 


from a calibration of the instrument. 


1.1-D Degrees of Freedom for Distillation Columns 

The available degrees of freedom for control of a distillation 
process is subject to varied opinion in the literature. The basis of 
the differences is the manner in which the various holdups within the 


column are treated. 
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Howard (65) states that the conventional steady state analysis 
technique is incorrect when applied to dynamical systems where the 
holdups of material may vary. The degrees of freedom for the unsteady 
State process may be derived in the usual manner by including the 
holdup in the same fashion as a normal stream. The holdup is analogous 
to a normal process stream in that it has the same variables except 
that a mass of material is involved instead of a flow rate. 

An alternate approach is given by Bertrand (66) and Murrill 
(67) who tabulate the variables and relationships for a distillation 
process. This approach considers only the steady state, the holdups 
being ignored. n+ 6 degrees of freedom will result from this approach, 
whereas the approach of Howard (65) shows n + 8 degrees of freedom. 
This difference, in the number of degrees of freedom, resulted because 
Murrill in his analysis assumed that the liquid levels in the condensor 
and reboiler were maintained constant by level control loops while 
Howard considered both levels as variables. Since the reboiler and 
condensor levels must be controlled the most logical approach seems to 
be to consider the column to have n + 6 degrees of freedom over and 
above specification of the reboiler and condensor levels. It should 
also be noted that if the approach of Howard is followed the top and 
bottom product flow rates could be controlled rather than the levels. 
If this is the case, the control of the holdup vessels would likely be 


over specified. 


1.2 PROBLEM FORMULATION 
As outlined previously the control of the composition of both 
the top and bottom product compositions from a distillation column 


presents special difficulties due to the interaction or dynamic 
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intercoupling that exists between the control loops. As a result of 
these interaction effects, the direct application of single loop con- 
trol theory will not result in satisfactory control of the column. 
There are basically two solutions to this problem of constructing a 
“two point" control system for a distillation process. One alterna- 
tive is the application of control systems as suggested by Rosenbrock 
(8) and Rijnscorp (5,6,7), the second alternative being the design of 
a multivariable control system for the process. 

It is the objective of this study to design and implement 
a noninteracting control system on a distillation column. The per- 
formance of the noninteracting system will be compared with a con- 
ventional two point beedback control system and the system suggested 
by Rijnsdorp (5,6,7). In order to design the noninteracting system, 
it will first be necessary to construct a simplified dynamic model of 


the process. 
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2. EXPERIMENTAL 


2.1 EXPERIMENTAL PROGRAM 

As outlined in Chapter 1 the objective of the present study 
is the comparison of systems for dual quality control of a distilla- 
tion column. This comparison was carried out by implementation of 
the control systems on a pilot scale distillation column. The syn- 
thesis of a noninteracting controller outlined in the previous chapter 
requires a knowledge of the dynamic behaviour of the process. In 
particular a linear transfer function model that relates the manipulated 
inputs to the process outputs is required. The first step of the 
experimental program was to perform dynamic tests on the process to 
establish the necessary dynamic relationships. After establishing a 
Suitable model for the process the control systems were implemented 


and their performance in rejecting disturbances compared. 


2.2 EXPERIMENTAL EQUIPMENT 


The experimental studies were performed on a pilot scale dis- 
tillation column interfaced with an IBM 1800 data acquisition and 
control computer. The distillation column employed in this study was 
an eight tray nine inch diameter column with a total condensor and a 
basket type reboiler. The column trays consisted of four 2-1/4 x 
1-7/8 inch bubble caps mounted in a square pattern. Solution holdups 
in the reboiler and on each tray were estimated by Svreck (68) at .25 
and .065 cubic feet respectively. For further details on the column 
construction the reader is referred to studies by Svreck (68) and 
Pacey (69). A compilation of operating procedures and equipment 


associated with the distillation column is on file with the Department 
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of Chemical and Petroleum Engineering (70). A binary mixture of methanol 
and water was the system used in the distillation column. The computer 
system provided the means for the rapid data acquisition and implementa- 
tion of the varied control systems necessary for the completion of this 
project. 

Since the dynamic models to be employed in the controller 
synthesis were assumed to be linear, it was necessary to select an 
operating point about which the dynamics of the column could be ade- 
quately represented by a linear mode]. The standard operating condi- 
tions about which the model was constructed are given in Table 2.1. 
Since the column was somewhat overdesigned, it was necessary to strike 
a compromise in the selection of the base operating condition. This 
compromise was necessary in order to maintain a reasonable composition 
profile throughout the column and at the same time obtain a measureable 
composition in the bottom product stream. The column also had a 
tendency to "flood" due to foam formation on the trays. A lower feed 
rate than one would expect from design considerations was necessary to 


eliminate the flooding condition. 


2.3 COMPOSITION MEASUREMENT 

In an earlier study on the same distillation column Svreck 
(68) obtained measurements of the compositions by measuring the 
capacitance of the solution. This method proved to be adequate for 
solutions of high methanol content but gave erratic results at low 
methanol concentrations. An alternate method, namely a gas chromato- 
graph, was selected instead of the capacitance method for measuring 


the bottom composition. 


asf fovdnos si nt beyoTams ed’ ot 2fatiom >tmeiyb eid soake 


ne Joalse oF yys22505i) 25W at isent! 80 oF bemUzz6 Saw atesiitays 
~sbe 3d bfvoo. mmutoo sni to astmsny dt is tdw ddods Iniog pntsevsa0 
-tbiion onrievsqo bvebist2 sty .Teboo sesoil # yd bstneesvqey ylateup 
F.8 sfdeT nt nevip sis botoiwenvo 2ew (Sabon sdf dortw tuods enott 
aiivt2 DY Wis2zs0s7 Zew JT ~benpresbysve Jeitwomos e5w nimyfoo aff some 
24dT  hotdibnoo portsvsgo 92ad Sid to norizelse ony AT setmoveyne. | 5 
noldreootes sfdsnoznss 6 ntatatem oF 1Sb90 nr yrseesosn saw 921 muqmo2 
afdsayuesgm 5 aletde oats atie2 alt ta brs nmilo> oft Juonquowis of for 
s bad oefe putifoy sd? smpayde tovbovg motiod Sit ni oF Faogmos 

best yawol A .2yend srt no notteuriot mao? of sub “bool?” o% gonebnag 
od ywe22599n 2aw 2norsetabrenca npftesb movi tosqxe bibow ano nett Stet 
-nottinnes pntboolt srt atantntts 


dooyve. mmloo norseltiteth omez ant fo inte yoi[yse ne i - 
ait) pntwuzean yd enotsieagnon of) Yo. ednamawuzsdn benitde (38) 

ot steupabs a oF boving’ wortson aay <norpator: ont Yo sonas faqs 
Wor 46 Sach aiden sve dud oat theintaat t enorautos 


18 


TABLE 2.1 


TYPICAL STEADY STATE OPERATING CONDITIONS 


Overhead product 
Reflux 

Bottom product 
Feed 

Steam 


Reflux 
Feed 
Steam 
Condensate 
Reboiler 
Plate ] 

2 


3 
4 
5 
6 
7 


8 
Condensor 


Flow 
(1bs./min.) 
1.18 
1.95 
la27 
2.45 
1.71 


Temperature 
Oy) 
1ol 
168. 
283 
Zeke 
209. 
AUIS Y 
194. 
hei 
Was 
164. 
156. 
T52. 
148. 
143. 


LOR Ole COl Oe a OY Oa © On ~ 


Composition 
(wt. % MEOH) 
96.0 
96.0 

=5 
46.5 


ae wot 
= (.prm\.2d!) 


0.88 . i ant 
O98 ) eet 
2, 3 KSI 

2, oP ah.§ 

- eT 
syutérsqnet 
re Ge 

\. lel 
0,.8af 
0.88 
e,\S8 


9.208 


19 


2.4 CONTROL LOOPS 

In order to insure that no outside disturbances influenced 
the dynamic tests and control] studies al] inputs to the column were 
monitored and/or controlled. A schematic diagram of the distillation 
column and control loops is shown in Figure 2.1. 

The flow rates and temperatures of the feed and reflux 
Streams were controlled as well] as the steam flow rate to the reboiler. 
To insure a constant pressure operation the condensor pressure was 
controlled by means of the cooling water flow rate and the differen- 
tial pressure across the tower was monitored. Liquid levels in the 
condensor and reboiler were controlled by manipulation of the top and 
bottom product flow rates respectively. Examination of the control 
system will show that eight variables have been specified which is 
consistent with the degrees of freedom analysis of Chapter 1. The two 
ljevels, the column pressure as well as the feed and reflux temperatures 
were controlled by analog controllers. Control of the top and bottom 
compositions was accomplished using the IBM 1800 computer. The required 
control actions (i.e. the reflux and steam flow rates) were cascaded to 


the setpoints of the flow controllers. 


2.5 COMPUTER INTERFACE 

The values of all variables described above as well] as the 
plate temperatures and temperatures of the flow streams were available 
for logging by the IBM 1800 computer system. Pneumatic signals were 
converted to analog signals suitable for the computer by air to current 
transducers. Control signals originating from the computer were con- 


verted to analog signals by a digita] to analog converter and transduced 
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to a current signal by current output stations (zero order hold devices) 
then further converted to pneumatic signals suitable for the control 
valve. A feature of the IBM 1800 computer was the IBM supplied direct 
digital control (DDC) program. With the exception of acquiring data 
for material and energy balance calculations, all communication with 

the process was done through this system. Each process control loop 
and data acquisition point is assigned a "loop record". The DDC pro- 
gram then allows the user to automatically read the analog input signals, 
calculate the contro] actions and output signals to the process at 
Specified intervals. The DDC program also had the capability of col- 
lecting and storing a Jarge amount of data over a specified time. 

This feature was used to advantage in acquiring data for construction 
of the process transfer functions. 

Due to the large amount of data required for the material and 
energy balance calculations a separate program was written which would 
read the desired analog input signals and convert these signals to 
engineering units for use in the material and energy balance calculations. 
This procedure was necessary in order to allow other users access to 


the facilities of DDC at the same time. 


2.6 GAS CHROMATOGRAPH 

The gas chromatograph used in the analysis of the bottoms 
composition was a modified version of the Beckman Series C Analyzer. 
The analyzer was equipped with a sampling valye of one microlitre 
capacity. Sample injection was initiated on command from the computer 
by means of an external contact output (ECO). Closure of the ECO 


activated a 115 volt relay which in turn activated a solenoid valve. 
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The solenoid valve controlled the operation of the sample valve by 
directing the application of air pressure to a diaphragm within the 
sample valve. Closing the ECO would thus result in movement of the 
1 microlitre sample volume from the carrier gas stream to the liquid 
sample stream. On opening the ECO a short time later the sample 
volume would return to the carrier gas stream and the liquid sample 
would be flushed into the chromatograph column. Liquid sample flow 
was maintained by a pump circulating solution from the reboiler 
through the chromatograph and back to the reboiler. A tabulation of 
the operating conditions for the gas chromatograph is given in Table 
(fae 

The 0 to 20 millivolt output from the chromatograph was 
read by the computer at one second intervals and stored by the DDC 
program for later analysis. The data acquisition phase lasted approxi- 
mately 150 seconds at which time the data was analyzed and a new cycle 
initiated. Details of the programs for the data analysis are given in 
Chapter 5. 
2.7 EXPERIMENTAL PROCEDURE 

During the performance of the dynamic tests the column was 
brought to the desired steady state using the local analog controllers. 
Once at the desired steady state the data acquisition phase was 
initiated at the computer. After approximately fifteen minutes of 
Steady state data was acquired the dynamic test was initiated by manual 
adjustment of the desired manipulated variable. The data acquisition 
was allowed to continue until the process once again reached a steady 
state, at which time the data was punched onto cards for off line 


analysis. 
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TAS EEC wc 


CHROMATOGRAPH OPERATING CONDITIONS 


Temperature 120s C 
Helium Pressure 30 psi 
Column Length eat 
Type Porapak Q 
Cycle Time 160 sec 
Retention Time Water 60 sec 


Methanol] 100 sec 


a 7 4 a - ny 
= pot, 
: 7 ; a 
- - 
: : 
3°0sT — BwtersqmeT 
: faq 0€ oo) mwwezert mut tot 7 
dt da . dyna nin fod 
- 0 weqe09 a | aay | 
Pas is 7 
ja2 04) : antT sfx 
‘982 0a ysIaW emlT nofsnatsh 


‘552 OOF «<  . Tomgatad . . 7 


24 


The procedure followed during the control studies was to 
start up the distillation column under analog control. Once the 
desired steady state was achieved the system was switched so that the 
analog reflux and steam flow controllers received their setpoints 
from the computer. The control algorithms could then be initiated 
and the process disturbed by changing the feed flow rate. 

The operation of the distillation column is illustrated on 
the McCabe-Thiele diagram in Appendix 3-A. The thermodynamic consis- 
tency of the temperature composition measurements are examined in 


Appendix 3-B. 
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3. DYNAMIC MODELS 


The greatest single difficulty in the application of modern 
contro] theory to a chemical process and in particular distillation 
columns is the construction of a suitable dynamic model for the system. 
It is well established that a distillation process can be described by 
a set of nonlinear material and energy balance equations. However, a 
model of such complexity is of little use in a process control system. 
The nature of the equations, especially the equilibrium relationship, 
make it difficult to obtain an adequate model by linearizing the des- 
cribing differential equations. In any case, if a model is constructed 
by linearizing a component material balance associated with each stage, 
it would then be necessary to reduce the order of the resulting state 
equation. There are well known techniques (71) for reducing the number 
of state variables in such an equation based on elimination of the 
States with the eigenvalues furthest from the origin. However, this 
approach results in a further approximation especially if the order of 
the state equation must be drastically reduced as would be the case in 
a large industrial column. In addition, the dominant eigenvalues may 
not necessarily be associated with the states of interest, i.e. the 
terminal compositions. Thus dynamic testing would appear to be a more 
reasonable approach to obtaining a simplified dynamic model than the 
theoretical approach. 

The nonlinearity of the equilibrium relationship would appear 
to make the strongest contribution to the nonlinear behaviour of a 
distillation column. This is especially true if the column is operated 


such that the end compositions are in the "pinch" of the equilibrium 
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curve. Composition excursions from the steady state towards the ends 
of the equilibrium curve would be severely limited whereas changes in 
the other direction would not be so limited. Thus in order to con- 
Struct a meaningful model from dynamic test data, it is necessary to 
restrict the magnitude of the input disturbance such that the process 
is not unduly effected by the inherent nonlinearities. This is valid 
if the model is to be used in a control system where the terminal 
compositions are under "tight" control. 

The procedure adopted in this study was to perform a set of 
dynamic tests on the distillation column in order to establish the 
relationships between the two end compositions and the inputs (reflux 
flow, steam flow, feed flow and feed composition). Since the object 
of these tests was to establish a process model about the specified 
operating point a pulse test was employed rather than a step input 
which would tend to yield information about the dynamics on excur- 
sions away from the operating point. A rectangular pulse was used 
for determining the model parameters. 

The process was assumed to be adequately represented by a 
combined first order time constant plus a time delay transfer function 
model relating the process outputs to each input. The model parameters 
were then determined by fitting the time domain response of the model 
to the recorded transient data. 

The "fit" was accomplished by employing Rosenbrock's direct 
search technique (72) using a least squares criterion. In order to 
simplify the search for the optimum parameters the process time delay 
was first calculated. This was accomplished by locating the point 


with the largest slope on the initial side of the response curve. The 
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time delay was then determined from the intercept of the line of maximum 
slope and the initial steady state value. 

The above procedure assumes a linear process behaviour about 
the operating point. The validity of this assumption and the use of the 
simplified model can be seen by comparing the model predictions with the 
experimental data. 

The resulting model parameters are given in Tables 3.1 - 3.4 
with the area under the associated input pulse. The values of the process 
gains are tabulated in (wt. %)/(lbs./min) except in the case of the feed 
composition runs where the units are (wt. %)/(wt. %). The time constants 
and time delays are tabulated in values of minutes. Further data con- 
cerning the dynamic tests (magnitude of the forcing pulse, and steady 
State operating conditions) are tabulated in Appendix 2A. 

The transient response of the process is compared with that 
predicted by the model for several tests in Figures 3.1 to 3.6. In all 
cases the process and model responses are in good agreement illustrat- 
ing that the simplified model adequately describes the behaviour of the 
system about the designated operating point. 

In order to arrive at a set of transfer functions for the pro- 
cess, the average model parameters were computed based on the tabulated 
parameters. The time delay associated with the bottom product composi- 
tion does not include the 2.5 minute analysis time of the gas chromato- 
graph, however, it does include the time delay in the sampling lines. 
Thus it is necessary to subtract an estimated .5 minute delay from the 
tabulated bottom product time delay to arrive at the true process trans- 
fer function model. The transfer function model of the process is given 


by equations 3.1 and 3.2. 


| alia | uaheiyrcarinet 
suods woiwsiad: e2anorg "1 b etmuzes yous oT ' 
atid to sen atid ban inna exdeial qpeneael : 
eg ditt enotiotherq ‘TSbom anit ent-ysanios Yd a. eles 
stab TetnouPreges 
b.€ — (6 eefdeT ni novip sxe 2xatameteg fobom paFafazen ait. ; 
22990719 sit to 2svlev aT .selug Jugat batetsazes’ ‘aid 8000 6ay3 aie tithe 
hast ant to S260 aii nt tqboxe (nfm\. cdl) \(2 dw) nt betelpdas avs. ants 
ztnetano> amrd ont  .(X-.dw)\GS dw) a8 atinu afd svadw emi nortteoqmea » he 
-noa Bigb qativut .29¢untw 7o etsy nt betetided sve ayshob  amtdubas. 
ybsese bis .sefuq onioret and To sbutinpam) cist otmenyb ont ‘entiriao 
AS xtbnogiA wh badatuded ste (eno httbnod oni drvaqe atede 
tend Atiw baysqnes 27 22e0e gots to senoqes) Sratenssd say i s 
[le nl .O1f of f.6 eenett nf etzst fenovee 107 febom oft ad eesaibeaq =~ 
-dert2ulti Insweops bope nt-ais aaznogesy (ahow big e2zsoovg ony 20289 
ond to Worversd silt cadtnaeab yfotsupshs fsbom betthiqmte sat Send gat 
tytoq piitts1s90 bsseaptaeb snd dy0ds madzy2 
-o1g Ss) vot 2fafiom? asters Yo tse 5 36 yin ot yabio nl ee 
bersi/tde? sit no bse5d boauqnes saw etesamesg habom spevavs ot »2290. 
-i20qm0d Sauboxq mostod sit Werw bay inane vsfab smd oT ey etsme6q 
~otsmoyHs 260 sit Fo omty zheylens afuitm ¢.S sig aby tant adiaaauall 
2anit pint tqnie. seit) nt \yetob ‘smks. ait aby font asoh 1? rewaword «rigor 
- ald mort Ysfob: stunTm 2. batsinties ns sosytidue of yvRzeaoan 27 
ta nual pve nrdnaisaaatat peptic _ 


28 


Ae|aq 
Ow | 


Ue} SU0) 
owl | 


NOILISOdWOD WOLLOS 


Aelaq 4U27SuU0) 
OWL | OWL | 


NOILISOdWOD dOL 
SLINS AG SISSL AS ids xX idaa 


late dled Et 


29 


ae oe Ss 
one SSL 
21 aq uUe1SU0) 
SUL | OWL | 


NOTLISOdWOD WOLLOS 


ie) Bes pe Cake 
Gos ye 
eke Lc Ode 
Bac vc nai, 
ea 7uUe}SU0D 
uLed OWL | OWL | 


NOTLISOdWOD dOL 


SLIASSY SLSAL ASIAd WVALS 


CocesavL 


ULey 


Pouy 
as[nd 


Lé 
Ge 
02 
02 


ON 
uny 


Bf 


a el 


30 


ead 3Ue4 Sud) 
OWL | OWL | 


NOILISOdWOD WOLLOE 


as 0°8 0°02 
7 b*Oll 2 Oil 
oe ee) 6°EL 


e[ ag Ue} SUO) 


ULed OWL | OWL 


NOTLISOdWOD dOL 


STN SAUSS 1 Saas Uild] Moa aGaa4 


Ge Sagas 


_ - : 7 : 
> a ri Pre 7 
| a Se SC 
—- 


a 7 7 
- 2 a a 


ss asueae erear 229 w 
a a " 


nhea mitt: 


_oyetae 
a 
FoF 
8 


3] 


5 : £°9 ea Coe Gigt= 36 

[3 Gh, O6L* €°6 cS? eal L8 22 

. fs = °9 Ca ev" ae LZ 

@ OL 8° OL L60° € 8 L 9¢ 70° LEl 02 

ea 7uUe}SU0D e [aq 7Ue}SU0D Pad ON 

OWL OWL | ULeD SUL | awh | ULed as[nd uny 
NOTLISOdWOD WOLLOG NOTLISOdWOD dOL 


SLINSAY SLSAL ASTAd NOILISOdWOD 0434 


ema avi 


ge 


0-0 417-5 35-0 52.5 
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FIGURE 3.1] 


PRODUCT COMPOSITION RESPONSE TO A PULSE IN REFLUX 
FLOW [RUN 22] 
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PLGUIRE sos 2. 


“PRODUCT COMPOSITION RESPONSE TO A PULSE IN 
STEAM FLOW [RUN 24] 
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PRODUCT COMPOSITION RESPONSE TO A PULSE IN 
FEED COMPOSITION [RUN 22] 
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FIGURE 3.4 


PRODUCT COMPOSITION RESPONSE TO A PULSE IN 
FEED FLOW [RUN 20] 
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FIGURE 3.5 


TOP PRODUCT COMPOSITION RESPONSE TO A PULSE IN REFLUX FLOW 
(TOP - RUN 21) AND STEAM FLOW (BOTTOM - RUN 25) 
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TOP PRODUCT COMPOSITION RESPONSE TO A PULSE IN FEED COMP. 
(TOP - RUN 21) AND FEED FLOW (BOTTOM - RUN 21) 
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The percent deviation of the model parameters from the mean 
value for each set of runs is recorded in Table 3.5. For the control 
inputs (reflux and steam flows) the agreement is good, the % deviation 


being less than 15% in all cases except one. 


: 12.8071 +5 -18.9e72°95 ; 
D 16.7s + J 21S eo 
= (ey) 
; 6.6e/°S -19.4e7£:98 ; 
B T0ds 1 Coy caren 
: 3.878: IS 22071 +78 - 
D TiC =| poe al 
| 
= (3.2) 
; 4.9 Sis ae? 28 : 
Bo Fes 1 ae ey F 


The larger deviations in the case of the top product composition gain 
to changes in reflux and steam flow rates can be attributed to the 
nonlinear behaviour of the column. From Tables 3.1] and 3.2 it can be 
observed that a somewhat larger gain results when the system is forced 
out of the "pinch" than occurs when the composition is forced into the 
endataros 

In the case of the disturbance inputs (feed flow and composi- 
tion), more data is required to arrive at an adequate model. Since the 
feed flow and feed composition data is superfluous with respect to the 
study this data will not be given further consideration. 


A comparison of the results of the pulse tests in Tables 3.] 
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and 3.2 show that similar model parameters resulted from pulses of 
different magnitudes and in opposite directions. Thus it would appear 
reasonable to assume a linear process behaviour with respect to the 
effects of reflux and steam flows on the end compositions over the 


small operating range considered. 


3.1 RELATIVE GAIN MATRIX 

The correct pairing of manipulated and controlled variables in 
an interacting control system may be accomplished through the use of a 
relative gain matrix (75). The object is to control a given process 
output with the manipulated input that will have the greatest influence 
On it. The relative gain matrix is determined from a "normalized" matrix 
of the open loop process gains. 

From equation 3.1 K is the matrix of the open loop gains. 
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An intermediate matrix L is then formed 

Reiko 
The relative gain matrix k! is then formulated by multiplying corresponding 
elements of K and L 

1.e. KG BS a) en pei 
For the case at hand the resulting matrix of relative gains between the 
manipulated variables (R and S) and the controlled process outputs (Xp 
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This illustrates that the overhead composition (Xp ) and bottoms 
composition (Xp ) should be controlled by the reflux (R) and steam (S) 


flow rates respectively. 
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4. DUAL PRODUCT QUALITY CONTROL SYSTEMS 


The most readily apparent scheme for accomplishing dual product 
quality control of a distillation process, shown schematically in Figure 
A. is control of the overhead composition by manipulation of reflux 
flow and bottom composition control by altering reboiler heating duty. 
However, it is evident on examination of the dynamic nature of the 
process (Equation 3.1) that a strong interaction will exist between the 
two control loops. Indeed the effect of the bottom product control on 
the overhead composition is larger than the direct effect of the mani- 
pulated variable, reflux rate, on the overhead composition. 

Thus any control system for the direct control of the overhead 
and bottom product composition must be based on a method that will com- 
pensate for, or reduce, the effects of the interaction between the 


overhead and bottom composition control loops. 


4.1 RATIO CONTROL 

A novel approach for reducing the effects of the interaction 
between the control loops has been suggested by Rijnsdorp et al 
(5,6,7). The control system that was suggested is illustrated in Figure 
4.2. In this scheme the overhead composition controller is cascaded to 
the setpoint of a ratio controller which controls the ratio of the over- 
head vapour rate to the reflux flow rate by manipulation of the reflux 
rate. A disturbance to the column which causes the overhead composition 
to change will cause a change in the desired ratio of the two flows, 
which for a constant overhead vapour rate will cause a corresponding 
change in the reflux flow. On the other hand the ratio controller will 


automatically compensate for a change in overhead vapour rate by a 
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FIGURE 4.1 


SCHEMATIC DIAGRAM OF THE CONVENTIONAL CONTROL SYSTEM 


4s 


FIGURE 4.2 


SCHEMATIC DIAGRAM OF THE RATIO CONTROL SYSTEM 
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corresponding change in the reflux flow rate. Since the effect of the 
overhead vapour rate on the top composition reflects the effect 

of the reboiler heating rate on the overhead composition, the top com- 
position loop has been rendered insensitive to control action in the 
bottom product loop. However, no account is taken of the effect of the 
reflux flow on the bottom product control loop. The effect of this 
omission on the performance of the control system remains to be 


evaluated. 


4.2 MULTIVARIABLE CONTROL 

A distillation column is a somewhat unique multivariable system 
due to the opposing actions of the two manipulated variables. Reflux 
flow has a positive effect on both end compositions; that is an increase 
in reflux flow will cause an increase in the more volatile components in 
both overhead and bottom compositions. Steam flow has the opposite 
effect as an increase in steam flow will cause a resultant decrease of 
the more volatile component in both end compositions. Consider a dis- 
turbance in feed flow or feed composition entering a distillation column 
where the top product composition is maintained constant by the reflux 
flow (top control loop) and the bottom product quality controlled by the 
steam flow to the reboiler (bottom control loop). Such a disturbance 
will have similar effects on both end compositions. If the disturbance 
is a step increase in feed flow for example, the top and bottom product 
compositions will increase. If the column is controlled as illustrated 
in Figure 4.1 the action of the controllers to maintain the end composi- 
tions at their required values would be to decrease the reflux rate and 
increase the steam rate. A decrease in reflux flow will cause a decrease 


in the bottom composition which is the desired action of the bottom control 
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loop. Similarly an increase in the steam flow will cause a resultant 
decrease in the overhead composition which is also the objective of 
the top control loop. However, the above analysis assumes that no 
oscillations will occur in the end compositions. If oscillations do 
occur the interaction between the control loops could oppose the effect 
of the manipulated variable and thus sustain the oscillations. There- 
fore, it would appear that it would be extremely difficult to tune the 
control loops if the interacting effects are ignored. The logical 
solution to the above proplem is to design a control system that will 
compensate for the interaction. This would serve to alleviate the 
difficulties of loop tuning and prevent one control loop from being 
disrupted by the action of the other. However, the above analysis has 
illustrated that the overall effect of the interaction, if controlled, 
is to enhance the action of the individual control loops. It may then 
be suggested that a more "optimal" method is one in which the interact- 
ing effects are used to enhance the performance of the control system. 

For the most part, the design techniques for multivariable 
control systems, reviewed in chapter 1, are concerned with servo control 
problem rather than the regulatory control problem. Only in a more re- 
cent study has Newell (76) applied optimal control theory to a regula- 
tory control system with unmeasured load disturbances. 

The approach for providing noninteracting control outlined 
by Zacklind (15) and Luyben (41), discussed in Chapter 1, appear 
to be the most applicable to the problem at hand. The basis of the 


method these 
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workers proposed is to compensate each control loop for the control 
actions taken in the other loops of the system. In essence feedforward 
control is applied to compensate for the effects of the interaction. 
The block diagram representation of the control system is 
shown in Figure 4.3. The process is assumed to be adequately repre- 


sented by 


x, (s) ; P,,(s) P.o(s) R(s) 
Xo (S) Po, (s) Poo(s) Sis) 


where the Pry are a combination of a first order transfer function and 
a dead time. The process model used in the synthesis of the control 
equation is given in equation 3.1. 

The design equations for the noninteracting compensator are 
constructed by the following analysis. The objectives of the compen- 
sator are to eliminate the effects of the top and bottom control 
actions on the overhead and bottom compositions respectively. 


Therefore 


Do, (s) Poo(s) + Po4(s) = 0 
and 


D,,(s) P+, (s) + P1o(s) = 0 


in general 


Dry(s) Pry (s) + Prgfs) =U 
Thus the compensators Diy may be constructed from the equation 


Diy (s) = - Py 3(s)/P 17 (s) 4.) 
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It is assumed in the above analysis that the dead time asso- 
ciated with the interacting function (P, 5) is greater than or equal to 
the dead time associated with the direct transfer function (Pry): 

This is necessary for the compensator Dry to be physically realizable. 
This assumption is usually valid as one would expect the dead time of 
the interacting transfer function to be greater than that associated 


with the direct transfer function. 
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9. IMPLEMENTATION 


The performance of the control systems described in Chapter 4 
were evaluated on the pilot scale distillation column. Implementation 
of these control systems was accomplished utilizing the IBM 1800 digital 
computer. The control computations were performed in programs executed 
in variable core on the 1800 at a specified periodic interval. The DDC 
program was employed as a convenient interface between the control 
algorithms and the process. The computed control actions to correct 
any error in the end compositions were inserted into the outputs of 
the specified DDC loop records. The outputs of these DDC loop records 
were cascaded to the setpoints of the analog flow controllers. All 
contro] computations were performed at 2.5 minute intervals as dictated 


by the cycle time of the gas chromatograph. 


5.) RATIO CONTROL SYSTEM 

The ratio control system studied was discussed in Chapter 4, 
in this system the bottom composition is controlled directly by a pro- 
portional plus integral (PI) controller manipulating the steam flow 
rate. The overhead composition on the other hand is controlled by cas- 
cading the output of the PI composition controller to the setpoint of a 
ratio controller. The ratio controller maintained the required ratio 
of overhead vapour rate to reflux flow rate. Since a direct measure- 
ment of the overhead vapour flow was not available, it was necessary to 
estimate the flow rate from a material balance about the condenser. If 
the level in the condenser is assumed to be constant, the vapour rate 
may be computed as the sum of the reflux and top product flows. The 


assumption of a constant level in the condenser is not unreasonable as 
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the level was under "tight" control during these studies. In addition, 
the period between control actions allowed sufficient time for the level 
to recover from any sudden changes in reflux flow made in the previous 
control interval. An average value of the top product flow rate was 
computed over the 15 seconds prior to the control action, using data 
accumulated, once per second, in a DDC ring buffer. Once an estimate 
of the overhead vapour rate was computed, the correct reflux flow could 


be computed on the basis of the specified ratio. 


5.2 NONINTERACTING CONTROL SYSTEM 

The block diagram representation of the noninteraction control 
system employed was shown in Figure 4.3. In this scheme the outputs of 
the proportional integral controllers are fed to the noninteracting 
compensators. The function of the compensators is to compute the neces- 
sary action to compensate for the effects of the single loop control 
actions on the other loop. The output to the process is the sum of a 
compensator output and the control action of the proportional plus inte- 
gral controller. To facilitate programming on the digital computer, it 
was desirable to express the required compensator (Equation 4.1) in a 
different equation format. Since a digital computer in a control system 
can be treated as a sampled data controller with a zero order hold 
device on the output, it is necessary to combine the process transfer 
functions of Equation 3.1 with the transfer function of a zero order 
hold in computing the sampled data form of the noninteracting compensa- 
tor. 


It follows that 
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In terms of the z transformed elements the compensator becomes 


Di y(Z) = -H Pry (Z)/H Pry (2) 
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where Wo is the compensator output and W. the input to the 
compensator. 


The above equation is restricted in that the delay time must 


be an integer number of sample periods. This is not desirable if a 
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long sample period is involved. To circumvent this restriction, the 
output was computed in two time steps. The time delay may be sub- 
divided into two intervals, the first interval representing the 
fractional sample period delay, the second representing the whole number 
of sample periods remaining in the time delay. The compensator output 
can then be computed in two time steps using the output after the first 
time step as an initial condition for the second time step. 

The numeric formulation of the noninteracting compensator 
and the ratio controller are outlined in Appendix 4. In addition the 
sensitivity of the estimated overhead vapour rate (by a material 
balance about the condenser) to changes in the condenser level is 


examined in Appendix 4. 
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6. PROGRAMS 


The control studies and associated data acquisition were 
performed by a number of programs executed on the IBM 1800 computer. 
For the purpose of description these programs may be divided into 
four categories (online programs, initialization and data retrieval 
programs, material balance programs, and auxilliary programs). A 


description of the auxilliary program category is given in Appendix 1. 


6.1 ONLINE PROGRAMS 
MAINLINE CHRTM 

The online programs performed the gas chromatograph data 
analysis and the control computations. These functions were controlled 
by the interrupt core load CHRTM, while the separate functions were per- 
formed by individual subroutines. 

On entry to CHRTM the ring buffer acquiring data from the gas 
chromatograph is deactivated and an ECO is closed to initiate a new 
sample injection to the GC. The core address of the ring buffer con- 
taining the GC data is located and the data loaded into a vector for 
analysis. Following computation of the composition as an area percent 
by subroutine GC the appropriate calibration is applied to convert the 
result to weight percent. The ECO is then opened to complete the sample 
injection phase and the ring buffer is activated. Subroutine CNTRL is 


called to implement the control action. 


SUBROUTINE GC 
Subroutine GC is called by CHRTM to analyze the accumulated 


data from the gas chromatograph. The start and end points of each 
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chromatograph peak are located by subroutine PEAK. The area under each 
peak is computed using Simpson's rule. If it is determined that the 
base line is reached at a peak end, base line corrected areas are com- 
puted for all peaks occurring since the last time an "on the base line" 
condition was detected. The base line corrected areas were computed 

by fitting a least squares line to adjacent sets of base line data. 

The superfluous area under the base line is then subtracted from the 
corresponding peak area. The corrected peak areas are assigned to 

the designated time bands and the area percent associated with each 


time band computed. 


SUBROUTINE PEAK 

Subroutine PEAK is called from subroutine GC to locate the 
start and end of each peak in the chromatogram. A peak start is 
recognized when the first derivative of the data becomes greater than 
a specified deadband. This event must be confirmed for a given number 
of successive data points (3) before a peak start is accepted. Similar 
logic is employed to locate the end of a peak; that is, the computed 
first derivative must become more positive than a specified deadband. 
This event must also be confirmed a designated number of times (3) 
before a firm decision is made. A test is made to determine if the 
base line is reached at the end of each peak. 

The derivatives used in peak detection are computed from the 
least squares formula given by Savitzky (64). The seven point formula 


used in this routine is 


y' = (3 Y(1#3) ¥ -2Y(it+2) 4+ Y(1*1) 
- ¥(I-1) - 2Y¥(1-2) -3Y(1I-3))/28H 
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where H is the time between successive data points. 


SUBROUTINE CNTRL 

Subroutine CNTRL was employed to interface the control algor- 
ithms with the DDC system and the calling program CHRTM. The subroutine 
was constructed in a manner which facilitated the use of the different 
control algorithms in the form of subroutines. The measurement of the 
bottom product composition is obtained from CHRTM and the overhead 
composition measurement from the DDC loop record. The composition 
measurements are computed as perturbation variables about the specified 
setpoints. The past history of control for the previous four sample 
periods is obtained from the disk file. 

The control algorithm is then applied. The controller out- 
puts are limit checked and converted to a percent value for use in 
the DDC loop record. The disk file is updated with the most recent 
measurements and controller outputs. 

As a proportional integral algorithm was used in all the 
studies, it was incorporated in subroutine CNTRL. The output of the 
PI controller was computed from 
tes KE eK } E. 

SUBROUTINE RCNTL 

The function of subroutine RCNTL is to perform the necessary 
control calculations associated with control of the overhead composi- 
tion by means of the overhead vapour rate to reflux ratio. Subroutine 
RCNTL is called from subroutine CNTRL. 

The desired ratio for control of the top product composition 


is obtained from the proportional integral control calculations performed 
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in CNTRL. An average value of the top product flow rate is computed 
from data obtained in a ring buffer. The overhead vapour rate is then 
computed as the sum of the reflux and top product flow rates. This new 
value of the overhead vapour rate is used to compute the new reflux 


flow rate necessary to yield the specified ratio. 


SUBROUTINE NICMP 

Subroutine NICMP was employed during the noninteracting control 
studies to compute the noninteracting compensation from equation 5.2. 
The output of the proportional integral controllers were fed to subroutine 
NICMP which calculated the compensating action based on these inputs. 
The control output to the process is then computed as the sum of a com- 
pensating term and a control action term. 

To circumvent the restriction of an integer sample period delay 
in equation 5.2 the compensation was computed in two time steps as out- 


lined in section 5.2. 


6.2 INITIATION AND DATA RETRIEVAL PROGRAMS 
MAINLINE CSTRT 
CSTRT was used to perform the following functions on request 
from the 1816 keyboard. 
1) Initiate repetitive execution of the gas chromatograph/control pro- 
gram CHRIM. 
2) Terminate repetitive execution of CNTRL. 
3) Initiate control action by setting a flag that informs CNTRL that 
the control action is to be implemented. The disk files are also 


initialized to start the control run. 
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4) Terminate control action. 


5) Change controller constants and/or model parameters. 


MAINLINE DUMP 
The nonprocess program DUMP was employed to retrieve the data 
from the disk files at the end of a control test. Measurements of the 
top and bottom compositions and controller outputs are dumped to the 
line printer or both the line printer and cards according to a 


specified option. 


6.3 MATERIAL AND ENERGY BALANCE PROGRAMS 
MAINLINE DASS 
DASS initiates the data acquisition for the material and energy 
balance computations. The current feed composition and bottom product 
composition are entered via the 1816 keyboard. The disk files are 
initialized and repetitive queuing of the data acquisition program DATAC 


at a specified interval is initiated. 


MAINLINE DATAC 

Program performs the data acquisition for the material and 
energy balance calculations. The program is queued a specified number 
of times at a given interval. A running total and sum of squares of 
the measurements are stored in the disk file. 

The data acquisition is accomplished by reading the multiplexer 
point directly and applying the appropriate conversion routines to the 
ADC reading. Three different calculation options are available, flow 
calculation, zero check, and thermocouple conversion (with cold junction 


compensation). These calculation routines parallel those in the DDC 


program. 
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The cold junction compensation and temperature conversion were 
accomplished following the example outlined in the IBM 1800 Functional 
Characteristics Manual (73). Thermocouple calibration data is stored 


as a series of straight line segments in the program. 


MAINLINE BALNC 
BALNC uses the data acquired by DATAC to complete the material 
and energy balance calculations. The average and variance of the 
accumulated data is computed. The flow rates are corrected for tempera- 
ture and composition. The material balance and a summary of the operat- 


ing conditions are printed on the line printer. 
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RESULTS 


The object of the control phase of the present study was to 
compare experimentally different methods for dual product quality 
control of a distillation column. The performances of the control 
systems discussed in Chapter 4, that is ordinary proportional integral 
control, overhead composition control by manipulation of the reflux to 
overhead vapour rate ratio and the noninteracting control system were 
studied on the pilot scale distillation column. The control loops 
were tuned for minimum integral of the absolute error (IAE). This 
enabled a comparison to be drawn between the ability of the different 
systems to compensate for a disturbance to the column. The control 
runs performed are summarized in Table 7.1. The steady state operating 


conditions for all control runs are documented in Appendix 2.1. 


7.1 SINGLE LOOP CONTROL 


The performances of the individual control loops that is 


proportional integral] control of only the top composition and proportional 
integral control] of only the bottom composition were studied. This 
allowed a comparison to be drawn between the performance of the single 
loop contro] with the behaviour of the system with both loops operative 
which results in interaction. The two control loops were tuned individ- 
ually for a minimum IAE response to a distuybance.. Figures /.li-and 7-2 
illustrate the performance of these two single loop control systems 
following disturbances in the feed flow rate. The performance of the 
single loop control systems even with the associated long sampling 
periods and dead time is quite adequate. The larger IAE term associated 


with the bottom composition control is due largely to the relative 
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FIGURE 7.2 CONTROL OF ONLY THE BOTTOM COMPOSITION [15% INCREASE 
IN FEED FLOW, RUN S-2] 
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instability of the bottom composition control loop. The rather "drastic" 
effect of the bottom product control on the overhead product quality is 
the most important observation to be drawn from these tests. This 
observation underlines the deleterious effect of the interaction on 


the controllability of the overhead product composition. 


7.2 PROPORTIONAL PLUS INTEGRAL CONTROL OF BOTH PRODUCT COMPOSITIONS 


From observations of the results in the previous section one 
would expect the interaction if neglected, to cause a serious deteriora- 
tion in the performance of the control system when an attempt is made to 
control both of the end compositions. In order to establish if indeed 
this was the case an attempt was made to control both of the product 
qualities using the controller settings determined for control of only 
One composition. The effect of attempting simultaneous control of both 
product compositions on the performance of the control system is 
immediately evident from Figures 7.3 and 7.4. The sustained oscillations 
that result are the predominant characteristic of the system's response 
to a disturbance. No IAE values are given for these tests as the system 
did not reach a steady state in a reasonable time. An IAE comparison 
would thus be meaningless. The disturbing effect of the interaction, 
steam rate, on the overhead composition is clearly evident from results 
given in Figures 7.3 and 7.4. As the overhead product composition, 
under the control of the reflux flow, begins to recover from the initial 
disturbance, a second disturbance effects the control loop causing a 
large deviation of the composition in the opposite direction to the 
initial disturbance. This second disturbance is the effect of the 
control of the bottom composition, i.e. steam rate on the top composition. 


The end result is a sustained oscillation. 
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FIGURE 7.3. PROPORTIONAL INTEGRAL CONTROL OF BOTH END COMPOSITIONS 
[15% INCREASE IN FEED FLOW, RUN PI-1 | 
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The performance of the system is clearly undesirable. On 
examination of the control behaviour of the loops it can be reasoned 
that a much "slower" control is necessary especially in the bottom 
composition control loop if the oscillations are to be reduced. On 
the basis of this reasoning an attempt was made with both control loops 
operating to tune the individual proportional integral controllers for 
a minimum IAE recovery from a disturbance. Figures 7.5, 7.6 and ta) 
illustrate the slow control that is necessary to achieve a satisfactory 
recovery from a disturbance. A complete tabulation of the controller 
constants employed in the control tests is given in Table 7.2. A 
substantial reduction in control action in the bottom loop was neces- 
Sary. The proportional and integral modes in the bottom control loop 
were 3.0 and 2.2 times less respectively than the corresponding con- 
Stants employed when only the bottom composition controller was active. 
This large reduction in bottom loop control action is necessary to 
prevent the overhead composition contro] loop from being overly dis- 
rupted by the interaction. Recal] that equation 3.1 shows that the 
overhead composition is more sensitive to steam rate than reflux rate. 
The reduction in controller constants is especially important in the 
case of the proportional mode and less important with regard to the 
integral term as sudden changes in the steam rate are less acceptable 
than a gradual change. 

In addition to reduced control action in the bottom control 
loop, it was also necessary to reduce the controller constants in the 
overhead composition controller to prevent a sustained oscillation. 

The associated proportional and integral constants were 1.7 and 2.2 


times less respectively than when only the overhead composition 
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PROPORTIONAL INTEGRAL CONTROL OF BOTH END COMPOSI- 
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PROPORTIONAL INTEGRAL CONTROL OF BOTH END COMPOSI- 
TIONS [15% DECREASE IN FEED FLOW, RUN PI-4] 


COMPOSITION (PERCT = ) 


COMPOSITION (PERCT +) 


7| 


OQ eS 55°0 62-5 110-0 
TIME (MINUTES) 


FIGURE 7.7 PROPORTIONAL INTEGRAL CONTROL OF BOTH END COMPOSI - 
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controller was active. 

A comparison of the IAE figures for the control runs made 
with only one control loop active against those obtained when both 
control loops were active clearly shows a deterioration in performance 
in the latter case. Indeed, the IAE approximately doubles when both 


control loops are active. 


7.3 RATIO CONTROL 

In this control system the overhead composition is controlled 
by means of the reflux to overhead vapour rate ratio. As discussed in 
section 4.] this ratio control of the overhead composition reduces the 
sensitivity of this composition to the interaction. It is at once 
evident from the results illustrated in Figures 7.8, 7.9, 7.10 and 7.11 
that the performance of this control system is a vast improvement over 
that attainable with conventional proportional integral control. It 
can be seen from these results that the overhead composition has been 
rendered insensitive to control action taken in the bottom loop. This 
observation can be drawn when Figures 7.8 - 7.11 are compared with 
Figures 7.3 - 7.7 where the overhead composition appeared to be affected 
by a second disturbance due to the interaction as outlined in section 
tee 

A comparison of the IAE's resulting from tests R-1 to R-4 
against those obtained with ordinary PI control in runs PI-3 to PI-5 
show that this ratio control system has succeeded in reducing the IAE 
associated with the top composition loop by a factor of about 2. In 
addition, an improvement is also noted in control of the bottom composi- 
tion. This improvement is due to the faster control that may be 


employed in the bottom control loop. This increase in control action 
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RATIO CONTROL OF THE TOP COMPOSITION [15% INCREASE 
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that may be taken in the bottom loop can be attributed to ratio control 
of the overhead composition. Care should be taken when comparing the 
performances of the different control systems on the basis of the IAE 
values for the bottom control loop as the composition analysis becomes 
unreliable below 0.3% weight percent methanol. 

A potential drawback of the ratio control system is that this 
system ignores the effects of control action in the top loop on the 
bottom composition. This does not appear to seriously effect the 
controllability of the column used in this study, but it should be 
noted that the reflux has a relatively small effect on the bottom product 
composition in this case. Also the stability of the bottom loop is such 
that tight control is impossible. This makes is difficult to evaluate 
the effect of the interaction on controllability in this loop. Thus 
it is difficult to evaluate the effect changes in reflux flow could 


have on the bottom control loop in other distillation control systems. 


7.4  NONINT 

The noninteracting contro] system was successful in achieving 
a significant improvement in control over the conventional proportional 
integral control system. The performance of this control system follow- 
ing a disturbance in feed flow is illustrated in Figures 7.12, 7.13, 
7.14 and 7.15. The noninteracting system reduced the IAE in the top 
composition control loop by a factor of 2 over that obtained with con- 
ventional control. In addition, there was a significant improvement in 
control of the bottom composition as the IAE in this loop was reduced 
from 29.3 to 16.7 in a typical case. 

In comparing the performance of the ratio control scheme with 


the noninteracting system, both systems result in a similar "goodness" 
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FIGURE 7.13. NONINTERACTING CONTROL [15% DECREASE IN FEED 
FLOW, RUN C-2] 
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FIGURE 7.14 NONINTERACTING CONTROL [15% DECREASE IN FEED 
FLOW, RUN C-3] 
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of control as illustrated by the IAE values for different tests. There 
is a marginal improvement in control in the bottom composition loop 
with noninteracting control over that obtained with the ratio control 
Scheme. 

A comparison of Figures 7.12 - 7.15 with 7.3 and 7.4 illustrates 
the effect of the noninteracting compensation. The overhead composition 
is no longer disrupted by control action taken in the bottom loop as 
was evident in Figures 7.3 and 7.4. The sustained oscillations evident 
in Figures 7.3 and 7.4 have been eliminated. In fact control of both 
end compositions together with noninteracting compensation yields as good 
a performance in the individual control loops as was attainable when only 
one end composition was controlled. 

The controller constants that yielded a minimum IAE response 
with noninteracting compensation were the same as those employed to 
obtain a minimum IAE response when two loops were controlled separately 
in runs S-1 and S-2. This suggests that noninteracting compensation 
allows independent operation of the two control loops. 

It is reasonable to suggest that noninteracting compensation 
would show an even greater improvement in the control behaviour of the 
column if the bottom composition control loop was more stable. The main 
objective of noninteracting compensation and the ratio control scheme 
is a reduction of the effect of control action taken in the bottom loop 
on the overhead composition. With a relatively unstable bottom product 
contro] loop, the amount that the contro] action in this loop can be 
increased with the use of an advanced contro! algorithm is limited. 

Thus the full benefit of the noninteracting and ratio control systems 


may not be realized in this case. 
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CONCLUSIONS 


A gas chromatograph has been demonstrated to be applicable to 

a feedback control system where the ratio of the sampling 
period to the process time constant was .17. 

As a design basis for a noninteracting control system, the 
dynamic behaviour of a binary distillation column (for excur- 
sions in product composition of up to 1. wt. percent about an 
operating point) may be adequately represented by a set of 
first order plus dead time transfer functions. Averaging the 
model parameters obtained from a series of on line pulse tests 
coupled with a fit of the transient responses in the time 
domain provided a satisfactory means of evaluating these 
transfer functions. 

The direct control of both the overhead and bottoms composition 
of a distillation column is not recommended unless a means of 
reducing the interaction is employed. It has been illustrated 
that the interaction, if ignored, can cause a serious deterior- 
ation in the performance of the control system. One must be 
prepared to accept either extremely slow control or an oscill- 
atory performance in such a system. In addition, as a direct 
result of the interaction between the control loops, it is 
unduly difficult to tune the individual controllers. 

Control of the overhead composition by manipulation of the 
reflux to overhead vapour ratio provides an effective means for 


reducing the interacting effect of the steam flow on the over- 
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head composition. The use of this ratio control system results 
in a significant improvement in the performance of the overhead 
composition control loop as compared to a conventional system. 
The control system designed on the basis of eliminating inter- 
actions yielded a much improved control of the overhead composi- 
tion and a somewhat improved performance in the bottom composi- 
tion loop over a conventional system. In addition, the problem 
of tuning the individual control loops has been greatly reduced, 
as the noninteracting compensation allows the loops to be tuned 
independently. 

In the case studied in this work, there is little difference 
between the performance of the ratio control system and the 
noninteracting system. The additional improvement in control 
with a noninteracting system may be offset by difficulties en- 
countered in model construction. On the other hand, if the 
qualities of multiple side streams are to be controlled by 
their respective drawoff rates a form of noninteracting compen- 
sation is necessary. In this respect a noninteracting control 
system will receive wider application if difficulties in model 


construction are not prohibitive. 


RECOMMENDATIONS FOR FUTURE WORK 


The construction of an adequate dynamic model is the single 
greatest difficulty in the application of modern control theory 


to a distillation process. Although much work has been done in 


=4giat cststabed to ae ef 
-Faannos baativsyo sit to Tovtnod bavorgint prdeinctig ion. 
=f2nqmon motsed. a3 nf aaiamatisd beverrqnt yerwamoze 6 bre ay 
noldoxg sit .nottibbs mf .usteye fsnatinevnna s yevo qoot nots 
sbsaunay ylseatp mad 26d bub t tovtno2 Taubivtbni ey pataut to 
benut sd ot 2q00! 344 ewolfs notisensqnos poisasysentaon oid 26 

vf tnabneqebat 
gonevatttbh sfitii 2t svorlt ,dtow etnd al betbuse gzeo oft nl 
ant bra: mat2y2 Tortnoa ofiay sft te Sonemioieq oft naswied 


fovtnos A? tnsmovorqmt fenoftibbs ett .meze¥e pnisosredatnon 


-n9 gartluarttib yd tsztte ad yam mot2y2 anftiosiedataon 6 dtw 
ait Tt ,dn6d vento oft nO .nottavytenos Fsbom nt beatnuos 

vd bal fovtnos od of 945 empsyie sbt2 siaitium to estitisup 
-naqnod pattoniwsininon to mot 6 estat Throws. svidzsqes7 afond 
fovtnen pnizssisininon « soaq2s) eri nl .yiseesoen ef nottee 
fsbom ni eats (ustttth 3F fotisatlage xobtw avisoan (fiw mete 
‘avittdtrions 3on 916 mottoussenoa 


——— a 


86 


this area a completely satisfactory method for determining a 
simplified dynamic model does not exist. As has been illustrated 
in the present study a set of transfer functions determined by 
on-line dynamic testing yields a satisfactory description of the 
dynamic behaviour of the process. However, dynamic testing is 
usually not desirable in an industrial process and is at best 
time consuming. Consequently, a method of simplifying the higher 
order dynamic models prevalent in the literature is required. 
There are well known techniques for reducing the order of a 

State variable system based on elimination of these states with 
the least significant eigenvalues. However, there is little 
evidence of an investigation of the accuracy of a model when the 
order of the system has been drastically reduced. Such an inves- 
tigation would be especially significant in the case of a distill- 
ation column where the eigenvalues are usually of the same order 
of magnitude. Thus a meaningful study would be an investigation 
of the representation of column dynamics by a reduced dynamic model 
of the process and subsequent study of its performance as a basis 
for implementing advanced control systems. 

The problem of formulating a dynamic model for noninteracting 
compensation can be largely overcome by using a simple gain com- 
pensation in the place of the dynamic compensator. The gains may 
be determined from simple step tests on the process or a proven 
steady model of the process. 

The ratio control system employed in this study was successful in 


reducing the effects of the interaction by manipulation of the 
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reflux to overhead vapour ratio. That is any change in the over- 
head vapour rate (which reflects the interacting effect on the 

top product) is automatically compensated for by a corresponding 
change in reflux flow. As was the case in this work, it is often 
impossible to measure the overhead vapour rate. This would suggest 
the use of an alternate system whereby the overhead composition is 
controlled by the top product drawoff rate, the condensor level 
being controlled by the reflux flow. Thus any change in the over- 
head vapour rate would be automatically compensated for by a 
corresponding change in reflux flow in order to maintain the con- 
densor level constant. A similar system may be employed in 
control of the bottom product composition whereby the composition 
is controlled by the product drawoff and the reboiler level is 
controlled by manipulation of the steam rate. This is the mater- 
ial balance control scheme proposed by Shinskey (75). 

Rosenbrock (8) and Davison (10) have described a method for two 
point control of a distillation column based on a model analysis 
technique as described in Chapter 1. It would appear worthwhile 
to investigate the performance of this system and compare the 
performance attained with that obtainable using other multivari- 
able control systems. 

As suggested in the review of the literature in Chapter 1, 

perhaps the control system with the most potential for dual 
quality control of a distillation process is one based on optimal 
control theory. Although it is difficult to formulate an 
"optimal" control policy for a regulatory system, investigation 


of this area would appear to be warranted. 
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NOMENCLATURE 
A coefficient matrix for the state equation 
aig element of the matrix A 
A analyzer 
B coefficient matrix for the state equation 
CG matrix of controller transfer functions 
Cij individual controller transfer function 
CC composition controller 
noninteracting compensator matrix 

ij noninteracting compensator element 

E error signal 

F distillation column feed flow rate 

FRC flow recorder controller 
G closed loop transfer function matrix 
Gi element of the closed loop transfer function matrix 
H(s) zero order hold transfer function 
Hq controller matrix 
K gain matrix 
Ki, individual gain element 
KS proportional gain 
Kr integral gain 
LC level controller 
n number of components in a stream 
lg process transfer function matrix 
Pay element of the process transfer function matrix 
R distillation column reflux flow 


R reference input vector 
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ratio flow controller 

distillation column steam flow 

sample period 

process time constant 

vector of control variables 

matrix of eigenvectors 

States yector 

distillation column bottom product composition 
distillation column overhead product composition 
distillation column feed composition 


matrix of eigenvalues 
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TRANSFER FUNCTION PROGRAM (TRSFCT) 


The function of the transfer function program was to fit the 
transient response data from a pulse test to a first order plus dead 
time transfer function model. 

The first function of the program was to smooth the transient 
data using an exponential filter. The transient data was then converted 
to perturbation values by subtracting a computed average steady state 
value. Since the process did not always return precisely to the 
initial steady state following the pulse test, it was necessary to force 
the pulse to close by insuring that the last 20 data points were indeed 
at the initial steady state. This procedure was necessary to avoid the 
large time constants that would result from any minor disturbance that 
would cause a steady state offset. 

To determine the parameters of the model the time delay was 
first determined by locating the point of maximum slope on the initial 
Side of the response curve. The time delay was then determined from 
the intercept of the line of maximum slope and the initial steady state 
value. The process gain and time constant were then determined by 
searching for the values of these parameters that minimized the sum of 
the squares of the errors between the actual response and that pre- 
dicted by model. The response of a first order transfer function pro- 


cess (=<) to a pulse input of width a and height b is given by 


mes a eee 
T 7 


Kb [1] - e - U(t - a)(J - e ] 


The search for the optimum process parameters was accomplished by employ- 


ing Rosenbrock's search technique (72). 


The program requires the following input data cards 
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1. NCASE. the number of cases to be run. 
2. For each case require card 1 + ICONT = 1 if the pulse input function 
is not part of the data input. | 
NPTS = number of points of transient data. 
ITYPE = 1 for a positive gain process. 
-1 for a negative gain process. 
NUMBR = number of vectors of data as input. 
Card 2 + NVCTR = 1 for each vector of data that is to be fit to a 
model. 
Card 3i% tittle 


Card 4 + PWDTH = pulse width in minutes. 


PHGHT = pulse height in process units. 
(IF ICONT = 1 the pulse width and height are calculated from 


the input data) 


IVCTR = vector number of the pulse input data (used only if 
ICONT = 1) 
NPSST = number of steady state data points, prior to the 


start of the pulse. 
NWIDE = pulse width as a number of data points. 
Card 5 + FLICST = filter constant. 
Cards 6 + + NPTS of transient data for each vector of data (preceded 


by 4 title cards which are ignored.) 


AUXILLARY TRANSFER FUNCTION PROGRAM (GCCON) 


This program is functionally similar to the transfer function 
program TRSFCT. The difference being that GCCON is employed to fit the 
transient data from the gas chromatograph to the first order plus dead 


time transfer function model. The data from the gas chromatograph was 
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necessarily acquired at a much slower rate than the other transient data. 

This necessitates several differences in the computer program to fit the 

transient data to the process model. 

1) A two point formula for the first derivative versus the five point 
formula used in the TRSFCT program was employed in calculating the 
deadtime. 

2) The last five points of transient data were forced to be at the 
initial steady state value in order to force the pulse to close. 

3) The transient data was not filtered. 


In al] other respects GCCONS is identical to TRSFCT. 


DATA SORT PROGRAM (GTDAT) 


The function of GTDAT was to sort the data acquired during the 
dynamic tests, convert the data to the appropriate units and punch the 
data to cards in a suitable format for the transfer function programs. 

Tne number of Joops, number of data points, time interval 
between data points and a flag (IGC) indicating that the gas chromato- 
graph analysis of the bottom product composition is to be outputted, 
for each of the data acquisition loops, the instrument range (A), 
offset (B) steady state temperature and composition are read from data 
cards. The steady state temperature and composition are used in con- 
junction with a flow loop to correct the flow rate for composition and 
temperature based on the steady state operating conditions. For other 
than flow loops the steady state temperature and composition should be 
Specified as zero. 

The data acquired by the DDC system is obtained from the data 
acquisition file by means of the system subroutine GBDAT. The data is 


converted from unscaled engineering units into actual engineering units 
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by the transformation. 


xX = A/32767 x Ix +.B 
where Ix is the unscaled value and x is the scaled value in engineering 
units. For flow loops the flow is corrected for the steady state 
composition and temperature. The data is then punched to cards. The 
gas chromatograph analysis is also outputted for the dynamic test if 


the IGC parameter is non zero. 
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APPENDIX 2 


A) STEADY STATE OPERATING CONDITIONS 2-1 
B) LOOP RECORDS 2-32 
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PULSE IN REFLUX FLOW(RUN 20) 


PULSE DATA 


PUESE HEIGHT 029 LBSe/MIN 


PULSE WIDTH 1467 MINUTES 


STEADY STATE OPERATING CONDITIONS 


TOP COMPOSITION 9565 WTe PERCTe 
BOTTOM COMPOSITION 051 WTe PERCTe 
FEED FLOW 2045 LBSe/MIN 
PEED, COMPOSITION 46066 WTe PERCTe 
REFLUX FLOW 1691 LBSe/MIN 
STEAM FLOW 2002 LBSe/MIN 
TOP PRODUCT FLOW 1618 LBSe/MIN 
BOTTOM PRODUCT FLOW 128 LBSe/MIN 


MATERIAL BALANCE ERROR OF CLOSURE 


OVERALL Oe6 PERCT.o 


METHANOL =O PERCTe 
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PULSE IN REFLUX 


PULSE DATA 


PULSE HEIGHT 


PULSE WIDTH 


STEADY STATE OPERATING 


TOP ‘COMPOSITION 
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FLOW(RUN 21) 
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CONDITIONS 
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PULSE IN REFLUX FLOW(RUN 22) 


PUESE DATA 


eeceecerecece 
PULSE HEIGHT 029 LBSe/MIN 


PULSE WIOTH 1363 MINUTES 


STEADY STATE OPERATING CONDITIONS 


TOP COMPOSITION 9600 WTe PERCTe 
BOTTOM COMPOSITION 055 WTe PERCTe 
FEED FLOW 2046 LBSe/MIN 
FEED COMPOSITION 4606 We PERCTe 
REFLUX FLOW 193 LBSe/MIN 
STEAM FLOW 2001 LBSe/MIN 
TOP PRODUCT FLOW lel7 LBSe/MIN 
BOTTOM PRODUCT FLOW 1230 LBSe/MIN 


MATERIAL BALANCE ERROK OF CLOSURE 
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PULSE IN REFLUX FLOW(RUN 23) 


PULSE DATA 


CCeeeeeccere 
PULSE HEIGHT 028 


PULSE WIDTH 8200 


STEADY STATE OPERATING CONDITIONS 


TOP COMPOSITION 9506 
BOTTOM COMPOSITION 050 
FEED FLOW 2046 
FEED COMPOSITION 4505 
REFLUX FLOW 1096 
STEAM FLOW ue rd 
TOP PRODUCT FLOW lei2 
BOTTOM PRODUCT FLOW 1030 


MATERIAL BALANCE ERROR OF CLOSURE 
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PULSE IN REFLUX FLOW(RUN 246) 


PULSE DATA 


eOoeeeseceve 
PULSE HEIGHT -e1l9 LBSe/MIN 


PULSE WIDTH 7070 MINUTES 


STEADY STATE OPERATING CONDITIONS 


@eeeeeevoeeeoevoe eee eevee eecesvseeaeveeseenoesee df 
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BOTTOM COMPOSITION 048 WTe PERCTe 
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REFLUX FLOW 1093 LBSe/MIN 
STEAM FLOW 2002 LBSe/MIN 
TOP PRODUCT FLOW lel7 LBSe/MIN 
BOTTOM PRODUCT FLOW le22 LBSe/MIN 


MATERIAL BALANCE ERROR OF CLOSURE 
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PULSE IN STEAM FLOW(RUN 20) 


PULSE DATA 


eeovoecesevece 
PULSE HEIGHT e200 


PUESE WIOTH 1007 


STEADY STATE OPERATING CONDITIONS 


@eeeaeetseoveeeeeeveeoseeoeeoeeeovoeeeseeeseeve eee 


TOP COMPOSITION 9506 
BOTTOM COMPOSITION 053 
PEED. FEOW 2045 
FEED? COMPOSTT ION 4606 
REFLUX FLOW 192 
STEAM FLOW 2002 
TOP PRODUCT FLOW 1220 
BOTTOM PRODUCT FLOW 1230 


MATERIAL BALANCE ERROR OF CLOSURE 
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PULSE IN STEAM FLOW(RUN 24) 


PULSE DATA 


PUES: GETor | ~e2l LBSe/MIN 


PULSE WIDTH veto MINUTES 


STEADY STATE OPERATING CONDITIONS 


e@eeeaoeeeoeoeeoeeeeoeseoeeooeseeseeeeegeoeeseone 


TOP COMPOST TION 9509 WTe PERCTe 
BOTTOM COMPOSITION 060 WTe PERCTe 
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REF GUX” PEOW 1696 LBSe/MIN 
STEAM FLOW 2001 LBSe/MIN 
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MATERIAL BALANCE ERROR OF CLOSURE 
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PULSE IN STEAM 


PULSE DATA 
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FLOW(RUN 25) 
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PULSE IN STEAM FLOW(RUN 27) 


PULSE DATA 


Ceooecececsce 
PULSE HEIGHT 


PULSE WIDTH 


STEADY STATE OPERATING CONDITIONS 


TOP COMPOSITION 


BOTTOM COMPOSITION 


FEED FLOW 
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STEAM FLOW 

TOP PRODUCT FLOW 


BOTTOM PRODUCT FLOW 
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PULSE IN FEED FLOW(RUN 20) 


PULSE DATA 


eCeeoeeeeevee 
PULSE HETLGHT 046 LBSe/MIN 


PULSE WIDTH 1562 MINUTES 


STEADY STATE OPERATING CONDITIONS 


TOP COMPOSITION hed | WTe PERCTe 
BOTTOM COMPOSITION 051 WTe PERCTe 
FEED FLOW 2045 LBSe/MIN 
FEED COMPOSITION 4606 WTe PERCTe 
REELUX ELOW 1094 LBSe/MIN 
STEAM FLOW 2002 LBSe/MIN 
TOP PRODUCT FLOW 1620 LBSe/MIN 
BOTTOM PRODUCT FLOW leat LBSe/MIN 


MATERIAL BALANCE ERROR OF CLOSURE 
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PUES IN] FEED “FEOWI RUN 24) 


PULSE DATA 


Ceoeseeeces 
PULSE HEEGHT 024 


PULSE WIDTH 155 


STEADY STATE OPERATING CONDITIONS 


@eeeooeaoeaeseveeeeeseeovoeeee eee eee eeeoeuvueeveeeoe 


TOP COMPOSITION 9509 
BOTTOM COMPOSITION 075 

PEED FLOW 2046 
FEED COMPOSTTION 4606 
REFLUX FLOW 1695 
STEAM FLOW 2000 
TOP PRODUCT FLOW 113 
BOTTOM PRODUCT FLOW 129 


MATERIAL BALANCE ERROR OF CLOSURE 
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PULSE IN FEED COMPOSITION(RUN 20) 


PUESE DATA 


eeoeeeresecs 
PUESE HEIGHT Bel WTe PERCTe 


PULSE WIDTH 1507 MINUTES 


SIEADY STATE OPERATING CONDITIONS 


TOP COMPOSITION 9545 We PERCT. 
BOTTOM COMPOSITION 052 We PERCTe 
FEED FEOW 2046 LBSe/MIN 
FEED COMPOSITION 4505 WTe PERCTe 
REFLUX FLOW 1294 LBSe/MIN 
STEAM FLOW 1099 LBSe/MIN 
(OPS rRODUCT FEOW 1.13 EBSe/MIN 
BOTTOM PRODUCT FLOW 1.28 LBSe/MIN 


MATERIAL BALANCE ERROR OF CLOSURE 
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PULSE IN PEED COMPOST TIONURUN 22) 


PULSE DATA 


eoeooeeeeeos 
BULSE HEIGHT “4069 WTe PERCTe 


PULSE WIDTH 1404 MINUTES 


STEADY STATE OPERATING CONDITIONS 


TOP COMPOSITION 9503 WTe PERCTe 
BOTTOM COMPOSITION 048 WTe PERCTe 
FEED FLOW 2046 LBSe/MIN 
FEED COMPOSITION 4565 Whe PERCT.« 
REFLUX FLOW 1094 LBSe/MIN 
STEAM FLOW 199 LBSe/MIN 
TOP PRODUCT FLOW 1.13 LBSe/MIN 
BOTTOM PRODUCT FLOW 1626 LBSe/MIN 


MATERIAL BALANCE ERROR OF CLOSURE 
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PULSE IN FEED COMPOSITTION(RUN 22) 


PULSE DATA 


@eeeeaeeoeeceeee 
PULSE HELGHT 509 WTe PERCTe 


PULSE WIDTH 1429 MINUTES 


STEADY STATE OPERATING CONDITIONS 


TOP COMPOSITION 9508 WTe PERCTe 
BOTTOM COMPOSITION 060 WTe PERCTe 
FEED FLOW 2045 LBS e/MIN 
FEED COMPOSITION 483 WTe PERCTe 
REFLUX FLOW 1696 LBSe/MIN 
STEAM FLOW 2000 LBS«/MIN 
TOP PRODUCT FLOW 1615 LBSe/MIN 
BOTTOM PRODUCT FLOW lel? LBSe/MIN 


MATERIAL BALANCE ERROR OF CLOSURE 
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PULSE IN FEED COMPOSITION(RUN 23) 


PULSE DATA 


eceecerevece 
PULSE HEIGHT —Tel We PERCTe 


PULSE WIDTH 1469 MINUTES 


STEADY” * STATE “OPERATING™® *CONDITIONS 


@eeeeeveseveeveeezs oeoeeeeeeseoeeoeaoeeeeaeenuee eo @ 


TOP COMPOSITION 9509 WTe PERCTe 
BOTTOM COMPOSITION 06 WTe PERCT.e 
PEED FLOW 2046 LBSe/MIN 
FEED COMPOSITION 483 WTe PERCTe 
REFLUX FLOW 196 LBSe/MIN 
STEAM FLOW 2000 LBSe/MIN 
TOP PRODUCT FLOW 1615 LBSe/MIN 
BOTTOM PRODUCT FLOW Lee LBSe/MIN 


MATERIAL BALANCE ERROR OF CLOSURE 
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CONTROL OF TOP LOOP ALONE(RUN $1) 


DISTURBANCE 


STEP IN FEED "FLOW 


STEADY STATE OPERATING 


+o 34 


CONDITIONS 


TOP COMPOSITION 


BOTTOM COMPOSITION 
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NONINTERACTING CONTROL (RUN C4) 


DISTURBANCE 


Sees Ne hob «PLOW 


STEADY STATE OPERATING 


CONDITIONS 


e@eeoeeeoeaeeveeeeevoveeeeeevneeeoeveeoevoeeeseegceeos 


TOP COMPOSI LION 


BOTTOM COMPOSITION 


PEED FLOW 

FEED, COMPOSI T.LON 
REPLUX FLOW 
STEAM FLOW 

TOP PRODUCT FLOW 


BOTTOM PRODUCT FLOW 


MATERIAL BALANCE ERRCR 


OF 


INITIAL 
96e1 


e660 


CCOSURE 


eecoeaeveeeeveeevneeveeeeaeeeoeeoeveoeoeeeoe eevee e sd 


OVERALL 


METHANOL 


—-003 


Oe9 


LBSe/MINe 


FINAL 
96el WTe PERCTe 


«60 WTe PERCTe 


280 LBSe/MINe 
4708 WTe PERCTe 
BER] LBSe/MINe 
ew unt LRSe/MINe 
1644 LBSe/MINe 


1641 LBSe/MINe 


pelt Cys PERCTe 


308 PERCTe 


a _ eonur gen #09 out en 
_ .- - a > _— - : 
. ry | an wa 
; ae - 
™’ |. We i a 
7 ia naes q re 
7 hod sponbe nee 
eMIM\ 9284 deen . — ona 183 if ay 
nd 7 ae _ 
: - _ “4 


: 3 ( i ss i. 
2O1TIGAO> an}TaRagO TATE ie 2 


CO eee 


JAALA SALT IMI an - 
& | Yas 
sTIASS «TW Lee Lea worTaoswa> sa io 


oT2A39. «Ti 08s Ode yorTy2o9Mo> oTTOR 


ies! 


EMA eBBD 8 08eS Tides . | 7 
; 8 dike he eee 
sT3#39 eTW | (Be Te Bete HOITI209M02 033% ; 
WIM\e@8i eeu abel -_ 
WIMN DORI) «6 FPaS Ses * 
' a de : == 
eAIM\ » QBS edel eSel : — WOJA T 
: A . = 
HIMN 2B bt ‘that wosa 


: - iti 


ET a 


- =) 
DES | Gig A. 
aa 


€ 


132 


DATA ACQUISITION LOOP RECORDS 


ID = 0170 BEED GhLOW 
6210 205C+04102 3340 8400407214. vUUu OCUUU Yuu 7EFF 
FPP 210 TFFF ©000 212106 


[ID = O17} REFLUX FLOW 
6U1U 205F+041U1 3340 8480408364 OOLU OULU VUUUU TFFF 
TERR 21 1lGelZGGGesoc00 2176 


ID = 0172 STEAM FLOW 
6010 205F+04100 3340 98480+05740 O0U0 OU00 ULvudO TFFF 
CEP “2ivOelgageecoooo —21L10 


ID = 0173 BOTTOM PRODUCT FLOW 
6010 205F+04104 3340 8480405072 ouLvU UU Uuuo TFFF 
7FFF 2110 7FFF 0000 2110 


ID = 0174 TOP PRODUCT FLOW 
6010 ZU5F+04103 3340 B4804+05382 0000 0000 @000 TREE 
7ePP 2110 T7FRF OO00 2110 


ID = 01975 TOP COMPOSITION 
6410 404C4+04106 O000 9100"01310409879 VJVOLU VOLO TFFF 
free 2a Freer OOOC Zale 


ID = 0176 BOTTOM COMPOSITION 
OZ1L0 4orc’ Cl77 GOOO0 8200420000 OOOO OOCO OOOO FFF 
Herr t2.20 FEF WOOO 2at0 


Fay 


FAT 


q44r 


qt 


33ar 


Jue 


Vid 


JULIU 


UYU 


UU 


VOU 


4 rr: 44 aa 
vuuy Luu arg vtoud G 
| olis, ue * 


a ae 

wosa xUUTEA ATT # sor 
uuu LUUG ABEBUSUBRE ORES fo,eoeneos ‘Urba: 
Gis guyu FAY Offs ee 


; his i 


WOR MAgT2 stro * a1 for 
GULU ULvuol OafeUebgeaR OFEE oatace3e0s ‘0408 ra 
OfIS OCov - <e OLIs 3 


ea’ Pay ae a Pa) 


:) Linea 7 gs A 
; Lag 


yuu ~LbLU SPOSUeuRaE iat SaoLso+szus 
‘7 if ofts heft ina Bi sf OLts 
Ln Le iets 


wO44 TOUGOHG MOTTOR —eyio = <1 a 
Rae 


wOJ> ‘TRudeR9 got 3=—ssi«éagO 
vuuu JJG SBEauUeUfas Usee . uke sper: a 
O1rs ooue Faar eT 


\ ' 


polmhaoewoa ser = ish LO 


Uy0o TBCU+OFETCs00L@ 0000 — sae dias fee 


T, me Waar 


a 
: pained 
ude eube a) 


Ql 
O02 


Ol 
O02 


Ol 
02 


Ol 
02 


Ol 
02 


Ol 
O2 


Ol 
02 


4010 
T7FFF 


0010 
7FFF 


4010 
TFFF 


7010 
TF FF 


7010 
TFFF 


7010 
7FFF 


7010 
TFFF 


133 


DATA ACQUISITION LOOP RECORDS 
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GAS CHROMATOGRAPH LOOP RECORDS 
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APPENDIX 3 


3-1) McCABE-THIELE DIAGRAM 


3-2) THERMODYNAMIC CONSISTENCY OF 
TEMPERATURE AND COMPOSITION DATA 
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McCABE-THIELE DIAGRAM 

The McCabe-Thiele diagram, Figure 3-1, is constructed 
employing vapour-liquid equilibrium data given by Perry (78). 
The operating lines for the two sections of the column can be 
written as: 


Section above the feed plate: 
R By 


> Ree r 


Yn ee eh 


substituting for R, D and XD from Table 2.1 yields: 
Yn =" 2625Xn +)" + 36.2 


Section below feed plate: 


R+F 


vite) ReaD 


X - BX 


m+ | B 


substituting for R, F, D, B, XB from Table 2.1 yields: 


Y = ATX =. 030 


m+ ] ner | 
(The temperatures of the feed and reflux streams were 
controlled to yield saturated liquids.) 


The McCabe-Thiele diagram illustrates the operation of 
the column in terms of theoretical stages. The extremities of 
the column clearly operate in the "pinched" region of the vapour- 
liquid equilibrium relationship. 
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100 


% Methanol in Vapour 


Wt. 


0 Wt. % Methanol in Liquid 100 


Figure 3-1 
McCabe-Thiele Diagram 
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THERMODYNAMIC CONSISTENCY OF THE TEMPERATURE COMPOSITION 
MEASUREMENTS 

When simultaneous composition and temperature measure- 
ments of a saturated liquid are available the measurements may 
be validated by comparison with published data. Temperature - 
composition data for a saturated solution of methanol and water 
(at 1 atmosphere total pressure) was obtained from Perry (78). 

In the case of this study the condensor effluent was 
subcooled leaving the reboiler as the only location where temp- 
erature and composition measurements of the saturated liquid 
were available. A comparison of the published data with the 
operating data shows a composition of 1.3 wt. % MEOH at the 
Operating temperature of the reboiler (209 .6°F) versus 0.5 wt. 
% MEOH measured in the reboiler effluent. (The operating 
pressure of the distillation column was -.1 inches of H,0.) 
Considering the purpose for which the composition measurements 


were used, a discrepancy of 0.8 wt. % is not unreasonable. 
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APPENDIX 4 


CONTROL ALGORITHMS 
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NONINTERACTING COMPENSATOR: 

The numeric forms of the noninteracting compensator diff- 
erence equations are derived by substituting the model parameters 
from equation 3.1 into equation 5.2. 

i Compensator Dl2: 
Wo(n) = 1.20 Wi (n-1.9) - 1.04 Wi (n-2.9) 
- .887 Wo(n-1) 

ii Compensator D21: 
Wo(n) = .437 Wi (n-4.1) - .368 Wi (n-5.1) 
- .796 Wo(n-1) 


RATIO CONTROL SYSTEM 

This control system controlled the overhead product compos- 
ition by manipulating the reflux to overhead vapour rate (V) ratio. 
The output of the proportional plus integral composition controller 
(0) was used to set the desired R/V ratio 


en R =@ 


from which the desired reflux rate was computed. 


The overhead vapour rate was calculated from a material 
balance about the condenser (a constant condenser level is assumed). 

ie. V = R+D 

The sensitivity of this calculation to a change in liquid 
level in the condensor can be evaluated from the change in condenser 
holdup per unit change in condenser level. 

Cross sectional area of condenser 


= .137 square feet. 
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Span of level instrument 
= 0.5 feet. 
Density of solution 
= 50 lbs./cubic feet (S.G. of 96% methanol water 
solution at 144°F = .80, ref. Perry (78)) 
Thus a change of 1% in the condenser level is equivalent 
to a change in liquid holdup of .034 lbs. This is equivalent to 1.1% 


of the normal overhead vapour rate. 
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